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We establish a relation between the two-point field strength correlator in QCD and the dual field propagator 
of the effective dual Abelian Higgs model describing the infrared behaviour of QCD. 



1. INFRARED ABELIAN DOMINANCE 

The QCD infrared dynamics appears to be well 
approximated by a (dual) Abelian Higgs model 
|TJ|^]. Via Abelian projection the QCD gluo- 
dynamics is reduced to Abelian fields, Abelian 
monopoles and charged matter fields degrees of 
freedom. Then, lattice simulations (mainly in 
SU(2) and in Maximal Abelian projection) show 
Abelian dominance and monopole dominance in 
the long range features of QCD. This means 
that nonperturbative dynamical quantities like 
the string tension a are reproduced (for more 
than 90%) simply by considering the Abelian 
or the monopole degrees of freedom. Moreover, 
monopoles are condensed in the confined phase. 
Electric fields and magnetic currents in the pres- 
ence of static sources are consistent with the 
dual Maxwell and Ginzburg-Landau equations. 
The penetration length A = 1/M (M = dual 
gluon mass) and the correlation length £ = l/M^ 
(At,/, — Higgs Mass) are measured. Finally, the ef- 
fective monopole Lagrangian is extracted from an 
inverse Montecarlo method and turns out to be 
equivalent to the Lagrangian of the dual Abelian 
Higgs model. Therefore, the above lattice investi- 
gations show that 1) QCD lattice gluodynamics is 
equivalent to the dual Abelian Higgs model: the 
Higgs particles are Abelian monopoles and con- 
dense in the confinement phase. 2) M ~ and 
therefore the QCD vacuum behaves as a dual su- 
perconductor on the border between type I and 
type II. 
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These are the information coming from the lat- 
tice. A still open problem is how to relate the 
Higgs and the dual gluon mass with the QCD 
parameters. In the following we show that an 
"Abelian effective dominance" shows up in the 
cumulant expansion of the quasi-static Wilson 
loop average in the nonperturbative region. This 
allows us to relate the non-local gluon condensate 
(i.e. the two-point field strength correlator) with 
the dual field propagator of an underlying dual 
Abelian Higgs model and therefore to suggest a 
connection between QCD and dual parameters. 

2. WILSON LOOP AND FIELD 
STRENGTH CORRELATORS 

In this section we report evidence of the Gaus- 
sian dominance in the cumulant expansion of 
the quasi-static Wilson loop average. The in- 
frared gluodynamics is contained in the Wilson 
loop which is an order parameter of confinement 

W(T) = P exp{ig <j) dz fl A tl (z)} and 1 is a loop 

made up by the two quark trajectories plus the 
endpoint Schwinger strings. In Euclidean space, 
using the non-Abelian Stokes theorem, the Wil- 
son loop average on the Yang-Mills action, ( ) = 

J VAexp{iSyyi[A]} , can be expanded as ||, 

(W(T)) = exp V ^ f dS(l) ■ ■ ■ dS(n) 

x (0(0, 1)F(1)0(1, 0) • • • F(n)<f>(n, 0)) cum . (1) 

where F(n) = i^ n i/„(u n ) is the field strength. 
S(T) denotes a surface with contour T and 



dS(n) = dS^ nl/n (u n ). The cumulants ( ) cum are 
defined as the connected part of the average and 
(0(O,l)F(l)0(l,O)) cum = (0(O,1)F(1)0(1,O)) = 
0. 0(0, u) = exp{ig J Q dz fl A fl (z)} is a Schwinger 
line. 
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((gF^(x)))r((gF X p(0))}r, with T a rectangular 
Wilson loop and the points x = (r, t) and be- 
longing to the first and second temporal quark 
line respectively. This object is required to com- 
pute the 1/m 2 corrections to the static quark- 
antiquark potential and has been evaluated on 
the lattice. The double bracket (( ))r stands for 
the average over gauge fields in presence of the 
Wilson loop W(T) (( )) r = ( W (V)) / (W (V)) . 
Taking into account Eq. ([l]), we obtain, 

(( 9 2 F^(x)F Xp (0)))r - ((gF^(x))) r ((gF Xp (0))) r = 
g {<f){Q,x)F liV {x)4>{x, 0)Fx p (0)} + 72^ ), (2) 

where 1Z contains all contributions coming from 
cumulants of order higher than two: 



Tlp, v \p(x;T) 
y> (ig) n 



n=3 {all perm. J a V- > 



j2 I n ds n» k ( u k) 

x(0(O, ui)F^ lUl (ui) ■ ■ -Fpuiui = x) 

x ■ • ■ F\ p (uj — 0) • ■ ■ Fp nVn {u n )(t>(u n ,Q)) cum . 



In the so-called Stochastic Vacuum Model |^] 
it is assumed that for large distances the bilo- 
cal cumulant (0(0, x)Fp, v {x)4>{x 1 0)F\ P (0)} is the 
dominant contribution to the r.h.s. in Eq. (||). 
This assumption has been (phenomenologically) 
successfully tested in potential theory and in soft 
high energy scattering processes. Eq. (||) is suit- 
able to obtain a first principles check of this as- 
sumption on the lattice. The adopted strategy is 
the following @. Within Eq. (§) 

Lorentz invari- 

ance is only violated by the contour T. There- 
fore, deviations of the l.h.s. of Eq. (||) from 
Lorentz invariance have to be attributed to the 
higher order cumulants 1Z on the r.h.s. We in- 
terpret the Lorentz invariant part of the data as 
(0(0, x)F flu (x)<j)(x 1 0)Fxp(0)) and parameterize it 
in terms of two form factors D and D\ 

( 5 2 0(O, a; )F^(x)0(x,O)F Ap (O)) = 



+ -g^-(x p Spx - xxS^p) 



D^x 2 ). 



(3) 



All the dynamics is contained in the form fac- 
tors. The function D is responsible for the area 
law and confinement. Both functions D and D 1 
are proportional to the gluon condensate G2 = 
(a s F 2 (0)}/n. 

In Ref. |H we have examined a set of lattice 
data for the l.h.s. of Eq.(]|). We have found 
that: 1) it exists a window (roughly x bigger 
than 0.2 -j- 0.3 fm) in which Lorentz invariance 
is restored in the continuum limit and therefore 
the r.h.s. of Eq. ^ is well approximated by 
the two-point field strength correlator; 2) in this 
window the functions D and D\ display an expo- 
nential decreasing behaviour ~ G2 exp{— |x|/T s } 
with a correlation length T g ~ 0.12 4- 0.19 fm. 
For the details we refer to Ref. j|. These results 
compare well with other lattice determinations of 
these quantities which have been obtained by use 
of the so-called cooling method || . However this 
is the first direct measurement without cooling 
and the first lattice evidence in favor of the Gaus- 
sian dominance in the Wilson loop average, 

(W(T)) ~ ^L^h^^^^ (4 ) 

In the long distance region (bigger than 0.2 fm) 
the Wilson loop average is well approximated by 
the two-point field strength correlator and there- 
fore the QCD infrared dynamics is "effectively 
Abelian" . Moreover, the nonperturbative regime 
is controlled by two parameters: the gluon con- 
densate (F 2 (0)) and the gluon correlation length 

T i 

given by 



g . The static quark-antiquark potential is then 



Vo(R) = \ I d 2 x {R-\ Xl \) D(x 2 )+^D 1 (x 2 ) 

Z J\x 1 \<R Z 

with d 2 x = dx\dx&, x 2 = x\ + x\. In the limit 
r/Tg 3> 1, Vo(r) — > or with the string tension 
given by 

° = \\ d2xD ^- ( 5 ) 



A non- vanishing D function leads to confinement. 

Quite naturally now a question arises. What is 
the relation between the "Abelian behaviour" of 
the Wilson loop average discussed above and the 
effective dual Abelian Higgs model suggested by 
the infrared lattice data mentioned in Sec. 1? 

3. DUAL EFFECTIVE DYNAMICS 

In Sec. 1 we have reported that the infrared dy- 
namics of QCD seems to be effectively described 
by a dual Abelian Higgs model with action 



S(C lt ,4>) = Jd 4 x ^G^(x)G^(x) 
+ l(D^)*(x)(D^)(x) + V((f>*(x)<P(x)) 



(6) 



where G^(x) = d t _ L C l/ {x) - d^C^x) + 
Gi^x), C M being the dual field and = 



dr 



\ a ^^(x 



y(r,a)) the 



lo Jo da dr _ 
Dirac field strength containing the quark cur- 
rents flowing around the Wilson loop. 4> is the 

A. 



Higgs field and V(tj>* 



\Y is the 



Higgs potential. The Higgs field is coupled to 
the gauge field C M via the covariant derivative 
D^cj) = (dfi + ieC^)4>. The parameters of the 
model are the dual coupling constant e = 2ir/g, 
the Higgs coupling constant A and the Higgs vac- 
uum expectation value <f>Q. The classical equa- 
tions of motion follow from the action (||). We 
emphasize that in this dual description the clas- 
sical approximation does apply. 

We have now two effective Abelian descrip- 
tions of the Wilson loop average: Eqs. (^)- 
(0) and Eq. (||). These are discussed in Ref. 
H. There the two-point field strength correla- 
tor (||) is related to the corresponding dual quan- 
tity Q ai \ P {x,y) = {5\ a 5 pl - 5 Xl 8 pa )8 A (x - y) - 
£n V \ P ^l3aajd^d^]C va (x,y), where JC pv is the dual 
field propagator, and therefore to the underlying 
dual superconductor mechanism. We summarize 
here the main results. 

3.1. London limit without quark sources 

Let us consider the pure vacuum without quark 
sources. In the case A — > oo we have <f> = 4>q and 



the Higgs field part of the action decouples. The 
dual field propagator K, pv = (C^C), satisfies the 
equation: 

{d 2 8 V f 1 -d v d p -e 2 cj)l5 vp )lC va (x, y) = -5 pa 6 4 (x-y) 
with solution (oo stands for London limit) 



K, 



lC°°(x,y) 



M 2 
M Ki(Mx) 



(2rrf 



with M — etfio the mass of the dual gluon. As a 
consequence we obtain 



G<j- f \p{x) = (5\ (7 5p~ ( - 6\ 7 Sp <J )D° (x 2 ) + - 



d 
dx x 

D?(x 2 ) 



with D°°(x 2 ) 

d 



M 3 Ki(Mx) 



D^(x 2 ) 



Air 2 x 

A—^K co {x 2 ). Therefore, assuming the London 
dar 

limit of a dual Abelian Higgs action as a model 
of infrared QCD, the predicted form of the two- 
point field strength correlator as in Eq. (|3j> , man- 
ifests a non-vanishing function D and an expo- 
nential fall off for the functions D, D\. Since 

/C°°(x 2 ) — > \ f * _i - e~ Mx + -.., the 
V ^Hoo2 (2^) 3 / 2 VM X 3/2 

dual gluon mass M plays the role of a correla- 
tion length. 

The logarithmic divergence of the string ten- 
sion calculated via (||) shows that the London 
limit is unphysical near quark sources where the 
Higgs field cannot be considered freezed to his 
vacuum value. In a physical situation we have to 
consider quarks and a Higgs field which vanishes 
on the Dirac string. 

3.2. Dual dynamics with quark sources 

Let us identify in the long range 

W(T) ~ ( e - s(c - 0) ), 

with the action S given by Eq. (||) and contain- 
ing quark sources via the Dirac string. Using the 
classical equations of motion we find 

S(C*A*) = \ ! dS ai {v) f dS Xp (u) 

1 JS{Y) JS(T) 



2W(x)) 2 



V(tf(x)) 



■k^I (R + I U- R l l - 1)) - It is also pos- 

in R 

sible to relate the Higgs condensate to the gluon 
condensate G2 — A</>q. 
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However, neglecting the dependence of the Higgs 
field, via the equations of motion, on the strings, 
taking into account the contribution coming from 
the Higgs part as a finite contribution to the 
string tension, we find that, also in the general 
case, Ga-yXp is equivalent to the QCD two-point 
non-local condensate (^) and in principle gives 
information on the validity of the decomposition 
(||) and on the existence and the behaviour of the 
functions D and D\. In the approximation in 
which the Higgs field depends only on the trans- 
verse coordinates, Q depends only on the parallel 
components. In this approximation we can eval- 
uate analytically the electric component of Q and 
find that its behaviour is controlled by a form 

factor D(x) = - — —Kq 



2%l 2 ' 




where I — y/SZ/SA, and S c — lim eC np (x±)/x± 
and Sa, = lim ed(x±)/xi are constants known 



^0 



in terms of the Ginzburg-Landau parameter 
A/e 2 . This suggests the identification of the corre- 
lation length T g with the dual quantity I. Notice 
that at variance with respect to the London limit 
result, here the correlation length is not simply 
given by the mass M of the dual gluon. 

The static potential can be calculated exactly 
without the use of an ultraviolet cut-off 



V (R) 



4tt R 2tt ./o 
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2w~P 



q 2 f H 

S c / dx^R-xx) 
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Higgs contr. 



= R^s c 



— > R — S c + R <th, 

R^oc ZTT 



(7) 



being the Higgs contribution to the static po- 
tential given by a linear term with string ten- 
sion (Th- Taking explicitly into account this 

contribution, the total string tension is a = 
2 

^—S c + (Th- In particular, for a superconductor 

Z7T 

on the border (A/e 2 = 1/2) we have Vq(R) = 



4. CONCLUSIONS 

In the assumption that the infrared behaviour 
of QCD is described by an effective Abelian Higgs 
model we have related the nonperturbative be- 
haviour of the two-point field strength correla- 
tor (g 2 F liU (x)(t>{x 1 y)Fxp{v)4>{y 1 x)) in QCD with 
the dual field propagator in the effective Abelian 
Higgs model of infrared QCD. In this way the ori- 
gin of the non-local gluon condensate is traced 
back to an underlying Meissner effect and the 
phenomenological relevance of the Gaussian ap- 
proximation on the Wilson loop is understood as 
following from the classical approximation in the 
dual theory of long distance QCD. In particu- 
lar the correlation length T g of QCD, which we 
know from direct lattice measurements, can be 
expressed completely in terms of the dual theory 
parameters. We have calculated analytically the 
static potential and the string tension. It turns 
out that the string tension is given by an integral 
over a function of the correlation length which can 
be identified with the non-local gluon condensate. 
There is no cut-off introduced in this calculation 
since it is not performed in the London limit . 
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